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ABSTRACT 

We show how to place constraints on cluster physics by stacking the weak lensing signals from 
multiple clusters found through the Sunyaev-Zeldovich (SZ) effect. For a survey that covers about 
200 sq. deg. both in SZ and weak lensing observations, the slope and amplitude of the mass vs. 
SZ luminosity relation can be measured with few percent error for clusters at z ~ 0.5. This can be 
used to constrain cluster physics, such as the nature of feedback. For example, we can distinguish 
a pre-heated model from a model with a decreased accretion rate at more than 5<r. The power to 
discriminate among different non-gravitational processes in the ICM becomes even stronger if we use 
the central Compton parameter which could allow one to distinguish between models with pre- 
heating, SN feedback and AGN feedback, for example, at more than 5<r. Measurement of these scaling 
relations as a function of redshift makes it possible to directly observe e.g., the evolution of the hot gas 
in clusters. With this approach the mass-Lsz relation can be calibrated and its uncertainties can be 
quantified, leading to a more robust determination of cosmological parameters from clusters surveys. 
The mass-Lsz relation calibrated in this way from a small area of the sky can be used to determine 
masses of SZ clusters from very large SZ-only surveys and is nicely complementary to other techniques 
proposed in the literature. 
Subject headings: 



1. INTRODUCTION 

Clusters of galaxies are powerful probes of cosmology. 
In particular, the cluster number density as a function 
of mass and redshift can provide unique constraints on 
the growth of cosmological structure and e.g. on dark 
energy. 

The next generation of Sunyaev-Zel'dovich (SZ) sur- 
veys, e.g., the Atacama Cosmology Telescope (ACT) 
and the South Pole Telescope 3 (SPT), will provide a cat- 
alog of thousands of clusters. The flux of the SZ is pro- 
portional to the integral of the den sity times the temper- 
ature of the hot gas in the cluster ijSunvaev fc Zeldovichl 
1980) and the flux limit of the catalogs will approxi- 
mately translate into a mass limit. However, the actual 
value of the cluster's mass for a given SZ flux depends 
on the details of the SZ-mass relation, which in turn is 
governed by cluster physics. 

Recent work has shown how this relation can 
be obtained by u si ng numerical simu la tions (e.g. 
IQn fc Rensonl ll2t)0l : Ida Silva"eTaTI {jS51; IMot.l et n ,l J 
1)2005!) : |Nagail (120051)'! or an alytical approxim ations 
Ce.g. IDqs Santos fe Dord llgOOl: iReid fe Spereel i2005l) : 
IRovchowdhurv et alJ <|2005MOstriker et all lj2005j) L Nu- 
merical models find that the SZ-mass relation is expected 
to be very tight, implying that clusters masses can be 
directly read out from the SZ observations once the SZ- 
mass relation has been calibrated. In addition, both nu- 
merical work and analytical approximations have shown 
that the integrated SZ luminosity of a cluster is rela- 
tively insensitive to input physics used in the models, 
such as pre-heating and energy injection by supernovae 
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or Active Galactic Nuclei (e.g. iMotl et all (|2005lk iNagail 
l|2l^ : lLl[^eTaLl EM) ). The same models, however, 
show a larger scatter with mass of the central SZ decre- 
ment parameter (i/q) and a much stronger dependence of 
2/o on the cluster physics. 

Ideally, we would like an independent and robust mea- 
surement of the cluster mass from observations. This 
would open the possibility to not just calibrate the SZ- 
mass relation directly from observations, but also to test 
intra-cluster medium (ICM) models, and therefore to di- 
rectly probe cluster physics. Arguably, weak gravita- 
tional lensing provides a direct way to measure the mass 
of clusters. The distortions of the background galaxies 
due to the cluster's gravitational field are sensitive to the 
mass along the line of sight, not just the intra-cluster gas. 
These distortions can be use d to reconstruct the clust er's 
projected mass density (e.g., [Kaiser fc Squires! EM3\ ). 

The weak lensing signal induced by galaxy clusters 
has been measured for a large num ber of systems (e.g., 
iClowe et all l)2005l) : iSheldon et all l)200H) and references 
therein). The amplitude of the lensing signal depends on 
the redshift distribution of the faint background galax- 
ies, which can be obtained from photometric observa- 
tions. Observations from the ground are limited by the 
atmosphere: both the maximum surface density of back- 
ground galaxies and the minimum error on the shear arc 
set by the seeing. While observations from space could 
circumvent this difficulty, there is no wide-field imager 
available in the near future to do this (ACT will cover a 
region of a few hundred sq. deg. and SPT 4000 sq. deg.). 
As both these factors set the minimum cluster mass that 
can be directly detected from weak lensing, a direct mass 
determination is possible only for fair l y mas sive clusters, 
w 1O 15 M (e.g.. lMarian fc Bernstein] <|2005[n . 

Here we circumvent this limitation by taking advantage 
of the wide multi- wavelength coverage provided by ACT 
and SPT with their planned extensive optical follow up 
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(Southern Cosmology Survey and Dark Energy Survey). 
By combining (stacking) the weak lensing signals from 
multiple clusters with roughly the same SZ luminosity, 
otherwise undetectable shear signal can be amplified, and 
thus an average mass determination can be achieved. 

The goal of this paper is to quantify the error in the 
calibration of the mass-SZ relation achievable by stacking 
weak lensing observations of SZ clusters. This also will 
allow one to investigate cluster physics since the scat- 
ter and amplitude in the SZ-mass relation can then be 
directly compared with predictions from numerical sim- 
ulations. In addition, a well calibrated mass— SZ relation 
can lead to a more robust determination of cosmological 
parameters from cluster surveys. 

Here, we will consider what can be learned by stack- 
ing clusters from a survey with the specifications of the 
two-year ACT SZ survey, assuming that the entire clean 
region (200 sq. deg.) is covered by optical weak lens- 
ing observations with a seeing of ~ 0.6 to 0.7 arcsec and 
20 source galaxies per square arcminute. The outline 
of the paper is as follows: in section 2 we compute the 
expected error in the mass estimator for the stacked clus- 
ters. In section 3 we explore how cluster physics can be 
constrained with the mass-SZ relation calibrated from 
the finding of section 2. We conclude in section 4 where 
we also outline the possible consequences of this approach 
for constraining the evolution of the hot gas in clusters 
and for obtaining a more robust determination of cosmo- 
logical parameters from clusters surveys. 

2. COMBINING LENSING SIGNALS FROM MULTIPLE 
CLUSTERS 

In ground-based experiments, the cluster weak lensing 
signal is smaller than the noise for masses less than 10 15 
Mq , due to the low number density of background galax- 
ies and the error with which their sheer can be measured. 
Here we illustrate how the weak lensing signal from mul- 
tiple lower-mass clusters can be combined to increase the 
signal-to-noise (S/N) ratio. In particular, we will con- 
sider measuring the average mass of several clusters in 
a given redshift bin and m ass range. We follow closely 
iMarian fc Bernstein! <)2005[) and define our average mass 
estimator as the weighted sum: 
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where d is the index of the stacked clusters, and 7t de- 
notes the average tangential shear of source galaxies in 
redshift bin s in an annulus with physical radius R k in the 
plane of lens d, and Wk s ,d are the weights, which should 
be chosen e.g., to maximiz e S/N. This is the same a s 
the mass estimator used bv IMarian fc Bernstein! ([2005) , 
but here it includes a weighted sum over different clus- 
ters, and the tangential shear is in physical rather than 
angular annuli. 
The variance of this mass estimator is 
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The weights Wk s ,d that maximize S/N are, 
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The term in square brackets is constant for all fc, s and 
d. Denoting this normalization constant by A, clearly 
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Like IMarian fc Bernstein! l|2005|) : 

iBartelmann fc Schneider! l)2001|) . we define 
Z(z s ,Zd) = ^jf^Q(z s — Zd) where 6 is the unit step 
function. The shear per unit mass for a hypothetical 
source at infinity is, 

- f \ lT(Rk,z Sl z d ) 

loo(0 k , z d ) = — (6) 

MZ(z s , z d ) 

The variance of jt from N galaxies in annulus fc and red- 
shift bin s is approximated by o 2 /N and tr 7 is assumed 
constant. We also assume a constant number density n g 
of galaxies per unit solid angle on the sky, distributed in 
redshift according to a Poisson distribution, 

-P(zs) = 7^z 2 s e^T . ( 7 ) 
zz 



The number of galaxies in an annulus of thickness dQ is 
27rn s smd d9, where 9 = R/Dd, so the number of galaxies 
in a physical annulus per unit redshift of a cluster at 

redshift Zd is 2im g V{z s ) sin (jj^j / D d dRdz s . 

Taking the continuous limit of the summation over fc, s 
in Eq. and R <C Dd, we obtain 

„2 
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J2 d Io°° dz s V(z 8 )Z*(z s , z d ) ^dR9jUR, z d ) 
where er^ is the contribution to the variance from a sin- 



;le len s at redshift Zd, as derived bv IMarian fc Bernsteinl 
( 2005). Thus, as intuitively expected, if one were to stack 
N identical clusters (clusters with the same mass at the 
same redshift) the density of background galaxies would 
increase by a factor N and thus the error on the mea- 
sured mass would decrease by a factor \/N. 

2.1. Cluster profiles 
We assume a NFW radial cluster profile (jNavarro et alJ 

= ( / ufe / ^2 - ( 10 ) 
[r/r s ){l + r/r s y 

We set p s and r s according to the clus ters mass and 
redshi f t, using the fitting form ulae from IBullock et alJ 
pOO"H) : iBrvan fc Normal |l99^ . 
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r s (M vir , 2) = 
C(M vir ,z) = 
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where G(z) is the linear growth rate and C is the concen- 
tration parameter, fi(z) is the ratio of the mean matter 
density to /9 cr (z), p C r(z) is the critical density of the uni- 
verse at redshift z, and a 2 (M) is the linear density field 
variance smoothed with a top-hat filter. The mass is 
related to the virial radius by 

47T 



M = M v 



-A(z)p CI (z)4 i 



3 

vir ' 



(17) 



The tangential shear is given by iWright fc Brainerdl 

(Hod3) : 

r s S c p cl - , \ 
1NFW = — g\x) 



4ttG D d D ds 
where x = :r. Thus 700 becomes 

3G ,D d (5 c , 



(18) 
(19) 

(20) 



A(D d )c 2 r 2 C 3 

Following iMarian fc Bernstein! (|2005f) . we take the up 
per limit of the second integral in Eq|S]to be Ru m = 2r vir 
obtaining 
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If we denote the numbers of clusters per unit redshift 
per unit mass per steradian as n c (M, z), and the frac- 
tion of the sky observed as / s ky, then the mass estimator 
variance becomes, 
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Lastly, we include th e lognormal scatter in the concen- 
tration parameter from lBullock e t al. (2001J), by defining 
the probability distribution, 



P(C',M,z) = 
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TABLE 1 

Mass bin sizes and number of clusters per bin for 
< z < 0.2 and 200 square degrees. 
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2.1.1. 



Computing the available cluster number 

In order to find out the mass error achievable for a 
realistic survey we need to compute how many clusters 
could be stacked for a given mass bin, redshift and sur- 
veyed sky area. We use the Sheth-Tormen mass function 
( Shct h fc Tormer ll999t) with E isenstein and Hu's analyt- 
ical fit l|ElserIsTern<^rrii¥l997l) for the linear mass power 
spectrum, and we normalize as to get 1000 clusters above 
2 x 10 14 M Q in 100 square degrees, as expected for an 
experiment with the specifications of ACT. 

To simplify our calculations we have found the follow- 
ing fitting formulae, 
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where A 1 = 1.19357, 61 
-0.02149738, and 



-0.42588341, c x 



da(R(M)) 
dR 
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where A 2 = 0.08502, b 2 = -0.4590054, c 2 = -0.0095691. 
The cosmological parameters we assume are il m = 
0.3, = 0.7 and h = 0.72 in a ACDM universe 
ijSpergel et al.l 12003^) . The above fitting formulae only 
apply for the cosmology chosen. 

2.1.2. Accounting for the effect of the large scale 
structure 
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M bin range (10 14 M Q ) Number of clusters 
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TABLE 2 

Mass bin sizes and number of clusters per bin for 
0.4 < z < 0.6 and 200 square degrees. 



M bin ranee (10 14 M Q ) Number of clusters 



0.50 - 0.62 


1241 


0.62 - 0.77 


873 


0.77 - 0.95 


602 


0.95 - 1.18 


405 


1.18 - 1.46 


265 


1.46 - 1.81 


168 


1.81 - 2.24 


103 


2.24 - 2.77 


61 


2.77 - 3.43 


31 


3.43 - 4.25 


18 


4.25 - 6.52 


13 


6.52 - 10.00 


2 



TABLE 3 

Mass bin sizes and number of clusters per bin for 
1.1 < z < 1.3 and 200 square degrees. 



The gravitational lensing signal is sensitive to all the 
mass along the line of sight. This is both a blessing and 
a curse: it makes lensing the most direct technique to 
measure clusters masses, but it also means that large- 
scale structure along the line of sight contributes to the 
lensing signal and affects the mass estimates. The effect 
of local large-scale structure (filaments connecting clus- 
ters and g roups) has been studied with numerical simula- 
tions ( e.g. ICenl l|1997j) : lMetzler et al1l|1999D : IWnite et all 
(2002)). This effect amounts to a bias in the mass deter- 
mination and an additional scatter. While there is not 



yet agreement on an estimate of the amplitude of the 
bias and th e scatt er induced by local structures (e.g., 
iClowe et alJ 1)2004(1 h we assume here that the bias can 
be accurately calibrated with numerical simulations and 
that the residual scatter is negligible compared to the 
weak lensing mass error. While this assumption is ex- 
pected to hold in most of the regime considered, it may 
break down at the high end of our mass range, thus more 
investigation is needed. This does not affect the results 
presented here significantly. In addition, the SZ-mass re- 
lation is expected to be affected less by these structures 
than the mass determination alone. In fact both weak- 
lcnsing mass and SZ signal are projections along the line 
of sight, the former of mass and the latter of temperature- 
weighted gas mass, and thus they are affected by these 
structures in approximately the same way. 

The effect of distant (not correlated with the cluster) 
large-scale structure is an important source of error in the 
mass dete rmination of a s ingle object, and it cannot be 
neglected. IHoekstral (|2001^) showed that these large-scale 
structu res introduce a n o ise in the clu s ter m a ss but not 
a bias . IHoekstral (|2001|) : iWhite et all l)2002fl : IHoekstral 
( 2003) quantify the effect and showed that it is roughly 
redshift independent for z > 0.1 and that it amounts 
to roughly doubling the error estimate obtaine d without 
accoun ting for the large-scale structure effect. IDodelsoiil 
l(2004j) showed how to somewhat correct for this effect, 
using the fact that the large-scale structure noise pro- 
duces a signal which is correlated over many pixels. This 
can significantly reduce the resulting uncertainties. Here 
we consider two cases: in the pessimistic case following 
IHoekstral (|2003|) we double the e rrors estimated a s in §2, 
in the optimistic case following [Dodelso rj |2004{) we in- 
crease the errors of §2 by 50%. In the text and figures 
we report the pessimistic case (in the optimistic case the 
errors will be reduced by a factor 1.15), but in the tables 
we report both. 

3. CONSTRAINING CLUSTER PHYSICS WITH 
MASS-SZ SCALING RELATIONS 

Here we will consider two quantities as a function of 
cluster mass and redshift: the total (integrated over the 
projected cluster area of a cluster) SZ luminosity, L sz , 
and the central Compton parameter, y . Both these 
quantities are expected to depend on clusters mass ap- 
proximately as power laws, with some intrinsic scatter. 
The details of the relation between mass and either L sz 
or yo depend on the physics of the Intra Cluster Medium 
(ICM). It is customary to parameterize these relation 
over a range of masses as power laws with two free pa- 
rameters (an amplitude and a slope of a linear fit in log- 
log space) and study how these parameters are expected 
to change for different cluster physics. Turning the argu- 
ment around by studying the effect of cluster physics on 
two parameters one is basically describing cluster physics 
with a two parameters model. Here we estimate how well 
these two parameters could be measured: using the er- 
rors calculated for the mass bins in the previous section, 
we estimate the error on the slope and amplitude of the 
mass-L S z scaling and the mass-yo scaling. To find out 
how this measurement can then be used to learn about 
cluster physics, we compare these findings with the slopes 
and a mplitudes for va r ious I CM physics models accord- 
ing to lReid fc Spergell l|2005[) . hereafter RS. RS examine 
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Fig. 1. — Mass error (<tm) vs. mass bin for three different redshift 
bins. The dotted lines are for < z < 0.2, the thin solid lines are for 
0.4 < z < 0.6, and the dashed lines are for 1.1 < z < 1.3. The thick 
solid lines denote constant S/N: the lower, S/N=10; the upper, 
S/N=l. The error &m has been doubled from the calculation based 
on the noise from galaxy shapes (sec. 2) in order to account for 
the effect of large-scale structure along the line of sight. The sky 
coverage is 200 sq. degrees, and the number density of galaxies is 
20 per square arcminute, distributed as a Poisson distribution in 
redshift with mean z = 0.99 (zo = 0.33). S/N will scale as the 
square-root of the sky coverage and as square-root of the number 
density of background galaxies. 



the mass-SZ scaling relations for several analytical and 
phenomenological models of the intra-cluster medium, 
which span the range of plausible models motivated by 
recent observations (see RS and references therein) . The 
approach presented here can of course be applied also to 
other analytical models and to predictions from numeri- 
cal work. We show which ICM models we can distinguish 
based on these two parameters. 

We examine the mass range from 5 x 10 13 M Q to 
10 15 M and consider 3 representative redshift intervals, 
< z < 0.2, 0.4 < z < 0.6, and 1.1 < z < 1.3. We 
will hereafter refer to these as redshift bins L (low), M 
(medium), and H (high), respectively. We choose the 
mass and redshift bins so that the error, ±om, is greater 
than the size of the bin, ±A^, and we have more than 
one cluster per bin. Note that the three redshift bins are 
disjoined, therefore for a fixed survey area stronger con- 
straints can be obtained by considering additional red- 
shift bins and/or combining together all redshifts. We 
choose the sky coverage to be 200 square degrees, and 
rig = 20 source galaxies per square arcminute from a 
Poisson redshift distribution with zq = 0.33 (mean red- 
shift 0.99). As the reported errors scale with and 

1/ \/ fsky , quadrupling the sky coverage or the galaxy 
density would halve the errors. Figure ^ shows the mass 
bins chosen and the error per bin for the three redshift 
bins and two lines of constant S/N. The mass bin sizes 
and the number of clusters per bin are reported in tables 

[3 and |31 In the figures we have included large-scale 
structure effects by doubling the errors (pessimistic case) . 

Of course, mass cannot be directly observed. Forth- 

4 Note that the error per bin depends on the bin size, bin sizes 
at a given redshift are approximately constant in log M except for 
the highest masses, but binsizes differ between different redshift 
bins. 



coming experiments (e.g., ACT, SPT) will observe the 
integrated SZ effect, L sz . In practice, therefore, clusters 
will b e binned by their L . z . From num erical simulations 
(e.g.. iMotl et~aTI pOOl : jNaeail pOO^ ) the M-L sz rela- 
tion is expected to be monotonic, L sz = A(M/M ) aM , 
with small (% level) scatter, enabling one to assign an 
average mass to a small L sz bin without introducing ad- 
ditional scatter 5 . One can of course test this assumption 
a posteriori by comparing the predicted mass error with 
the intrinsic r.m.s. in a given bin for bins with high 
enough S/N. 

Here we assume a fiducial relationship between mass 
and L sz , with a M = 1-66 and A = 6.18 x 10~ 6 Mpc 2 
at M Q = 1.39 x 1O 14 M0, corresponding to a self similar 
model, without non-gravitational physics (i.e. heating 
and cooling) where e.g. the gas temperature is given 
solely by the dark matter virial temperature. 

Using this fiducial model, we thus transform the mass 
bins into bins in L sz . The results are shown in figures 
and0]for mass bins L, M and H respectively. We then 
compare these scaling relations a nd their errors with th e 
theoretical models worked out bv lReid fe; Spergell {2005), 
hereafter RS. 

Our fiducial model corresponds to their self-similar 
phenomenological model. RS use a different equation 
for the concentration parameter C(M,z), but since we 
average over a range of C, it amounts to a small effect 
which we neglect here. Their models focus on the en- 
tropy parameter K = P/p 5 ^ 3 , where P is pressure and 
p is the density of the gas, so the gas entropy is propor- 
tional to IrtK. The gas entropy profile is parameterized 
by a normalization (-ftT ma x) and a radial profile with a 
power law index (si), and optionally a core with smaller 
power law index. Cooling and/or preheating affect K max 
and the profile. Preheating the gas causes it to start 
with a greater "preheating" entropy K p h, which serves 
to help offset the radiative cooling. Models of radiative 
cooling can reproduce important features of trends in X- 
ray clusters. However, including only radiative cooling 
cooling produces too much cooled gas, so some feedback 
is necessary. Other models deal with varying the accre- 
tion pressure. In the RS analytical models, the pressure 
is computed from the ram pressure of the infalling gas. 
The accretion is assumed to be smooth and spherical, 
while in reality is lumpy and stochastic. Deviations from 
these assumptions would change the accretion pressure. 
RS also examine models whose entropy profiles are de- 
rived by assuming a simple polytropic equation of state, 
P oc p 1 , with pre-shock gas density p\ and post-shock 
gas density pi (the gas is shocked as it accretes onto the 
cluster, when it crosses the sound barrier). 

Our fiducial model corresponds to RS self-similar 
model with an entropy profile K cx r , a baryon frac- 
tion of 0.13, an entropy normalization -Kniax = -ftTioo- 
The two models with extreme values of ctu are the RS 
analytical model with the lowest amount of preheating, 
-Kpreheat = 10 30 erg cm 2 g -5 / 3 (and the largest fraction of 
cooled gas, / coo i = 0.42), which we will call model 1, and 

5 The reader may note here that the scatter in the M-L sz relation 
can introduce something similar to Malmquist bias as lower mass 
clusters are more numerous than higher mass ones. Once the width 
of the scatter in the Wl-L sz relation is known this effect can be taken 
into account. For our purposes, given the small scatter found in 
simulations, we neglect this (small) effect. 
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Cluster model 




A (kpc 2 ) 


Fiducial (self-similar) 


1.66 


6.18 


Model 1 


1.77 


4.76 


Model 2 


1.64 


5.74 


Model 3 


1.71 


7.78 



TABLE 4 

The best-fit slopes a M and amplitudes A for RS models 

PLOTTED IN OUR FIGURES, WHERE L sz = A{M/M ) a " AND 

M = lO 14 M0/h. Model 1 is the RS model with preheated 

ENTROPY PARAMETER K = 10 30 ERG CM 2 G" 5 / 3 , MODEL 2 IS THE 
RS MODEL WITH DECREASED ACCRETION RATE aj accr = 1/3.5, AND 
MODEL 3 IS THE RS MODEL WITH INCREASED ACCRETION RATE 
^accr — 3.5. 



the phenomenological model with an accretion pressure 
decreased by a factor of 3.5 from self-similar spherical col- 
lapse, model 2. The models with extreme values of ampli- 
tude A are model 1 and the phenomenological model with 
an accretion pressure increased by a factor of 3.5 from 
self-similar spherical collapse, model 3. Table 0] shows 
the best-fit slope, au> and amplitude, A, for the M-L sz 
relation for these models. We have included these three 
models in the figures, where the dotted-dashed line de- 
notes the low preheating value (model 1), the dashed line 
denotes the lower accretion pressure (model 2), and the 
dotted line denotes the higher accretion pressure (model 
3). Performing a joint analysis using both the slope and 
and the amplitude reveals that we can distinguish mod- 
els 1 and 3 at the > 3tr level in the low redshift bin, 
and at better than 5 a in redshift bin M. For the survey 
parameters and redshift bins chosen, we cannot distin- 
guish model 2 very well in any redshift bin (see table |5J). 
However, if we combine measurements from all redshifts 
or if we increased the size of the survey to 1000 sq. deg. 
and, in this case, consider only redshift bin M, we could 
distinguish it at the 3cr level. 

We expect that some systematic effects may affect 
much more one parameter than the other. If this is the 
case, one may want to consider the constraints on only 
one parameter (marginalized over the other one) . There- 
fore, next we look at what can be learned using just slope 
or just amplitude information. 

For the slope, we find that in the pessimistic case the 
error on the slope is a aM = 0.095 for bin L, a au = 0.052 
for bin M and a aM = 0.24 for bin H. 

We can thus distinguish model 1 from the fiducial to 
about 2a in redshift bin M, but cannot distinguish much 
else in this case from slope information alone (see table 

EJi. 

Since the error in the slope is so small for redshift bin 
2, due to the good signal to noise (and therefore large 
number of mass bins allowed), we can also distinguish 
among closer RS models in this case. For example, in 
the pessimistic case, we can distinguish between the two 
phenomenological models with varying accretion (by fac- 
tors of 3.5 and 1/3.5) to more at about 2a, and we can 
distinguish between the two extre me slopes of t he RS 
preheating models to more than la. iNagail (|2005|) found 
A a a/ ~ 0.1 between a purely adiabatic cluster model and 
one that included cooling and star formation. Clearly in 
redshift bin M for a sky coverage > 200 square degrees, 
we can begin to distinguish differences of this order and 




Fig. 2. — Left: Recovered mass vs. input mass and Lsz for 
redshift bin L: < z < 0.2. The horizontal lines denote the mass 
bin and the vertical lines denote the error on the average mass in 
that b in. Right: Using the self-similar model from Rcid & Spcrgcl 
12005) as a fiducial model (solid line), we have transformed each 
mass bin into a bin in L sz , and again plotted the bin width and 
error as crosses. The dotted-dashed line shows model 1, the dashes 
line shows model 2, and the dotted line shows model 3. Note that 
the slope in this plot is 1/cxm . 

i a Input Moss (M s ) , c- 
10 10 D 
1.25T ' — " I 
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Fig. 3. — Left: Recovered mass vs. input mass and Lsz for 
redshift bin M: 0.4 < z < 0.6. The horizontal lines denote the 
mass bin and the vertical lines denote the error on the average 
mass in that bin. Right: Using the RS self-similar model as a 
fiducial model, we have transformed each mass bin into a bin in 
L S z, and again plotted the bin width and error as crosses. The 
dotted-dashed line shows model 1, the dashes line shows model 2, 
and the dotted line shows model 3. Note that the slope in this plot 
is l/a M -l/a™ UCIAL . 



thus discriminate between these models. 

Both analytical and numerical work (e.g , JNagai 
(2005 ): IReid fc Spergell ll2005T) : IMotl et all {2005); 
iMcCart hybrid] l)2003fl ) show that different cluster 
physics modifies the amplitude of the M — L sz scal- 
ing much more than the slope. Nagai finds a differ- 
ence of about 30% in amplitude between his purely adi- 
abatic simulation and one with cooling and star forma- 
tion. RS models differ by up to 90% in amplitude at 
M = 10 14 M Q /h. 

Thus we are very interested in how well we can con- 
strain the amplitude of this relation. To do this, the weak 
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Cluster model 




# cr from fiducial, pessimistic 


# <t from fiducial optimistic 






z bin L 


2 bin M 2 bin H 


z bin L 2 bin M 


z bin H 


K ph = 10 30 (Model 1) 


3 


>5 1 


3.5 >5 


1.5 


Waccr = 0.28 (Model 2) 


0.5 


1.3 <1 


1 2 


<1 


Waccr = 3.5 (Model 3) 


3.5 


>5 1 


4 >5 


1.3 



TABLE 5 

The number of a FROM the fiducial to which the RS models plotted in our figures can be distinguished, using a joint 

ANALYSIS WITH BOTH SLOPE AND AMPLITUDE, IN THE THREE REDSHIFT BINS AND IN BOTH THE OPTIMISTIC AND PESSIMISTIC CASES. 
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Fig. 4. — Left: Recovered mass vs. input mass and L sz for 
redshift bin H: 1.1 < z < 1.3. The horizontal lines denote the mass 
bin and the vertical lines denote the error on the average mass in 
that bin rescaled by the fiducial. Right: Using the RS self-similar 
model as a fiducial model, we have transformed each mass bin into 
a bin in L sz , and again plotted the bin width and error as crosses. 
The dotted-dashed line shows model 1, the dashes line shows model 
2, and the dotted line shows model 3. Note that the slope in this 
plot is l/a M - l/af] DUCIAL . 



Fig. 5. — Cluster physics has a bigger effect on the the scal- 
ing of mass with yo. This plot shows assorted lines from RS 
for the mass-j/o scaling, with amplitudes normalized to agree at 



M 



10 Mq/1i. The crosses show the mass bins and errors for 



the fiducial model, the RS self-similar phenomenological model, for 
redshift bin M. The gray regio ns show the intrinsic scatter in yo 
according to Mo tl et alJ 120051) for the two most extreme slopes, 
models 4 (dotted-dashed) and 5 (dotted). 



lensing determination of clusters masses need to be accu- 
rately calibrated to yield the clusters virial masses. Here 
first we assume that this can be achieved with a negligi- 
ble residual bias, although this is indeed an assumption 
whose validity remains to be tested. Then we show how 
relaxing this assumption (i.e. allowing a residual system- 
atic error) degrades the results. 

Assuming gaussian errors and marginalizing over the 
slope, we find the la constraints on the amplitude A of 
the M - L sz relation to be 5.67 - 6.55 x 1(T 6 Mpc 2 for 
redshift bin L, 5.85 - 6.34 x 1CT 6 Mpc 2 for redshift bin 
M, and 5.3 - 7.0 x 1(T 6 Mpc 2 for redshift bin H. 

Using amplitude alone, we can distinguish models 1 
and 3 from the fiducial at better than 3cr, especially in 
the M redshift bin (see tabled). In the pessimistic case, 
we can distinguish all of the preheating models (from 
K = 5 x 10 33 to K = 10 30 ) from the fiducial to more 
than 2a in redshift bins L and M. 

Up to now we have assumed that the calibration of the 
mass determination has negligible residual systematic er- 
ror. A systematic uncertainty in the mass determination 
of x% would propagate into ~ 1.66 x x% systematic in 
L sz to be added to the uncertainties in the amplitudes 
reported above. Thus, for example, x — 5 would halve 
(reduce by a factor 2.5) the number of as reported for 



bin M. 

As poin t ed ou t by, e.g., iReid fc Spereell l|2005|) : 
iMotl et ail l)2005|h the scaling of the central Compton 
parameter, yo, with mass is far more sensitive to the 
cluster model than the integrated L sz . The effect of 
cluster physi cs on the L RZ — U p rela tion has been stud- 
ied by e.g. IMcCarthv et al.l (|2003). Here instead we 
investigate what additional information can be learned 
from measuring yo m addition to Lsz and weak-lensing 
mass determination. Even with a single mod el, however, 
there is a large scatter in the M — yo relation l)Motl et alJ 
12005(1 . This partly counter-balances the advantage that 
the slopes for different models are farther apart. 

As argued above, clusters could be binned in L sz , and 
the average mass per bin computed from the stacked clus- 
ters using a fiducial model. Then one could measure the 
yo for each cluster in the bin and take the average. This 
will have an additional error associated because of the 
intrinsic scatter in the M — yo relation. 

Taking the RS self-similar model as our fiducial, with 
a slope 0.84 and an amplitude 4.92, we plot the expected 
bin sizes and mass errors for redshift bin M as crosses 
on the fiducial line (see figure EQ). W e show the ex- 
pected scatter from IMotl et all l)2005l) . who computed 
the scatter from simulations of about 100 clusters, as a 
shaded gray region around the extreme models. Note 




yo 



Fig. 6. — Assorted models from RS for the mass-yo scaling, 
with amplitudes not normalized. The crosses show the bins and 
errors for the fiducial model, the RS self-similar phenomenological 
model, for redshift bi n M. The gray regi ons show the intrinsic scat- 
ter in yo according toJMotl ct al. (200Ef) for the two most extreme 
amplitudes. The dotted-dashed lines show the extreme preheat- 
ing models, the dotted lines show the extreme polytropic models, 
the dashed lines show the extreme accretion pressure models, and 
the triple-dotted dashed lines show the different values for sj, the 
exponent of the radius in the radial profile of K. 



that for a survey of 200 square degrees, one expects to 
detect about 2000 clusters, so the effect of this scatter 
on the determination of slope and amplitude of the scal- 
ing relation should be reduced. The two RS models 
with extreme slopes for this scaling are the analytical 
model with the largest preheating (-Kprehcat = 5 x 10 33 
erg cm 2 g -5 / 3 ), model 4, and the polytropic model with 
7 = 1.2, pi = 1.06pi, model 5. We can distinguish model 
4's slope from the fiducial model at more than the 5cr 
level for all redshift bins, and can distinguish model 5 
quite well in the lower two redshift bins (see tabic [2J - 

Figure shows the lines for various RS models without 
normalized amplitudes. The models with extreme ampli- 
tudes are model 4 (the preheating model) and model 6, 
the model with 6 si = 1.5. We can distinguish both mod- 
els' amplitudes from the fiducial model at 3> 5 sigma in 
all redshift bins and we can start constraining the pa- 
rameters describing cluster physics, e.g., the amount of 
entropy injection. 

Table |21 shows a few RS models, their amplitudes for 
the yo — M scaling, and their a from the fiducial. Notice 
we can easily distinguish models with amplitudes that 
vary by about 20%. Physics that causes such a differ- 
ence in amplitude includes changing the entropy profile 
normalization by a factor of 1.5, doubling the radius at 
which the entropy maxes out, changing the accretion rate 
by a factor of 3.5, or increasing the exponent S\ of the 
radial entropy profile from 1.1 to 1.5. 

As the slope for the yo relation is 0.84 (range 0.71 - 1.5), 
an additional x% systematic error on the mass determi- 
nation would propagate into an additional ~ 0.8 x x% 
error on the amplitude of the scaling relation. Even al- 
lowing for a 10% systematic mass error we find that, for 
redshift bin M, we can distinguish all of the RS preheat- 

6 Recall that si parameterizes the entropy profile, K oc r si and 
that si = 1.1 for the fiducial model. 



ing models, polytropic models, si = 0.7, and Si = 1.5, to 
more than 5tr. We can distinguish the models with vary- 
ing accretion pressure (models 2 and 3) to more than 2a. 
The range of amplitudes A y are 7.0 x 10~ 6 to 9.4 x 10~ 5 
(fiducial 4.9 x 10~ 5 ). 

Since the models with maximum and minimum slope 
and amplitude are different for yo than for L sz , combin- 
ing the two would allow for further distinctions among 
models of cluster physics. 

In the cases where the two parameters fit yields a 
S/N > 5 in principle there could be enough S/N to add 
parameters to the fit, or in other words to investigate if 
e.g., deviations from a power law in the scaling relation 
can yield additional information about cluster physics. 

Here we have concentrated on a survey with the prop- 
erties and sky coverage of ACT and its optical follow up. 
Naturally a wider SZ survey with weak lensing follow up 
and/or a weak lensing follow up from space will improve 
these constraints. A space telescope can yield n g = 100 
galaxies per square arcminutes. While for low and inter- 
mediate redshift this will improve constraints reported 
here only by y/E, as the additional galaxies are likely to 
have a distribution with a higher-redshift tail, the con- 
straints at the higher redshift will show a better improve- 
ment. This may open up the possibility to study redshift 
evolution of the ICM. On the other hand a ground-based 
survey covering 1000 square degrees will also yield con- 
straint better by a factor y5, providing better discrimi- 
native power between models. Finally a combination of 
ground-based SZ observations over 1000 square degrees 
with space-based weak lensing follow up would yield im- 
proved constraints by a factor 5. This combination is 
not too dissimilar from the planned SNAP observations 7 
if ACT could scan in SZ the same region of the sky. 

4. DISCUSSION AND CONCLUSIONS 

We have shown how, stacking the weak lensing signals 
from multiple clusters with roughly the same SZ lumi- 
nosity, the SZ luminosity-mass (L sz — M) and SZ central 
decrement-mass [y a — M) relations can be measured from 
forthcoming SZ surveys with extensive optical follow up. 
For example, from a survey of 200 square degrees such as 
ACT, we find we should be able to determine the scaling 
slope to ±0.05 — 0.2 and amplitude to 5-20% depending 
on redshift bin. Observations indicate that the entropy of 
the gas in clusters cannot be explained by gravitational 
collapse alone, but the cause of the increased entropy 
remains to be determined. Different non-gravitational 
processes affect both the slope and amplitude of these 
relations, and in particular the y — M relation is most 
sensitive to cluster physics. Used in conjunction with re- 
sults from simulations, the expected errors on the mea- 
surements of the slope and amplitude of scaling relations 
imply that we can start to discriminate among different 
non-gravitational processes affecting the ICM. 

We have shown that, for a survey of ~ 200 sq. deg., the 
available S/N enables one to measure these scaling rela- 
tions as a function of redshift: we can distinguish slopes 
and amplitudes of the M — L sz relation and M — yo rela- 
tion for different redshift bins with width Az = ±0.1, 
thus opening up the possibility to constrain how the 

7 snap.lbl.gov 
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Cluster model 


OLM 


# <*a M 


from fiducial, pessimistic 


# Ca M 


from fiducial optimistic 






2 bin L 


z bin M z bin H 


2 bin L 


2 bin M 2 bin H 


K ph = 10 30 (Model 1) 


1.77 


1.1 


2 0.4 


1.4 


2.6 0.6 


oj accr = 0.28 (Model 2) 


1.64 


0.5 


1 0.2 


0.7 


1.3 0.3 



TA BLE 6 

TWO EXTREME SLOPE MODELS WE CONSIDER, DESCRIBED BY ReID & SPERGEj] 120051) . THEIR BEST-FIT SLOPES Ct M (WHERE 
Lsz = A(M/M ) aM ), AND THE NUMBER OF a BETWEEN THEIR SLOPES AND THE FIDUCIAL SLOPE IN THE THREE REDSHIFT BINS AND IN BOTH 

THE OPTIMISTIC AND PESSIMISTIC CASES. 



Cluster model A (kpc 2 ) 



# a A from fiducial, pessimistic # a a from fiducial optimistic 







2 bin L 


2 bin M 


2 bin H 


2 bin L 


2 bin M 


2 bin H 


K ph = 10 30 (Model 1) 


4.76 


3 


6 


1.8 


4.5 


8 


2.3 


<^accr = 3.5 (Model 3) 


7.78 


3 


5 


1.6 


3.9 


7 


2 



TAB LE 7 

TWO EXTREME AMPLITUDE MODELS WE CONSIDER, DESCRIBED BY REID & SPERGEL (2005), THEIR BEST-FIT AMPLITUDES A (WHERE 
L 3 z = A(M/M ) aM ), AND THE NUMBER OF crS BETWEEN THEIR AMPLITUDES AND THE FIDUCIAL AMPLITUDE IN THE THREE REDSHIFT BINS 

AND IN BOTH THE OPTIMISTIC AND PESSIMISTIC CASES. 



mass-SZ scaling evolves with redshift. A mass measure- 
ment from weak lensing and a SZ measurement in differ- 
ent redshift bins, can constrain the evolution of hot gas 
as a function of redshift, which in turn would enable one 
to constrain feedback evolution. Inspection of Fig. Inland 
0] provides a quantitative measurement of how physics in 
clusters, e.g. feedback, can be constrained as a function 
of redshift. The errors reported here will scale with the 
survey parameters as l/-^// s ky and Thus, using 

this method, wider surveys can improve constraints on 
cluster physics just as well as deeper surveys. For a space 
telescope with n g = 100 galaxies per square arcminutc 
and 1000 square degrees, our results would improve by a 
factor of 5. 

A calibration of the SZ-Mass relation with known un- 
certainties also has consequences for the determination 
of cosmological parameters. The number density of clus- 
ters within a given mass range or above a given mass 
is extremely sensitive to cosmology. This sensitivity, in 
principle, can lead to tight constraints on cosmological 
parameters such as f2 m , dark energy properties and neu- 
trino mass. However, even relatively small systematic 
errors in clusters mass estimates can lead to system- 
atic errors on the recovered cosmologi cal parameters that 
are larger than t he statistical ones iHolder et alJl200ll 
iBattve fc Wellerll2003l iFrancis et alJl2005|) . Thispoten- 
tially compromises measurements of parameters like e.g. 
the dark energy equation of state. The technique out- 
lined in this paper to calibrate the mass-SZ relation may 
help reduce the amplitude of these systematic errors. 
The weak lensing mass determination used here is sen- 
sitive to the mass along the line of sight, not just the 
gas. It is true that there can be biases introduced by 
the large-scale structure local to the clusters, but we ar- 
gue that this should not be considered a fundamental 
limitation as this effect in principle can be quantified by 
means of numerical N-body simulations. If this can be 



achieved, our estimates above indicate that the resid- 
ual uncertainty in the calibration of the amplitude of 
the mass-SZ relation is 5% for a survey with n g = 20 
galaxies per square arcminute and sky coverage of 200 
square degrees. This uncertainty would propagate into 
a systematic error in the determination of cosmological 
parameters, that is now reduced to be not larger than 
the statistical unc ertainty achievable f rom a survey of 
the same size (see IFrancis et all l)2005j) L Alternatively, 
The M-Lsz relation calibrated in this way from a small 
area of the sky can be used to determine masses of SZ 
clusters from very large SZ-only surveys, or the M — yo 
relation could be used to constrain cluster physics and 
this information can then be used to model the M — L sz 
relation. This approach is nicely complementary to other 
techniques proposed to calibrate the mass-observable re- 
lations such as the self-calibration using the cluster power 
spectrum of Maiu mdar fc Mohrl l|2004ft . While only posi- 
tions and redshifts of clusters are needed to measure the 
power spectrum and thus to self-calibrate the SZ-Mass 
relation, it requires observations of large areas of the sky 
as square and contiguous as possible. The calibration 
of the M-SZ relation with weak lensing has more strin- 
gent follow up requirements (imaging with good seeing 
in addition to clusters redshifts), but has less stringent 
geometry requirements. 
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z bin L 


z bin M z bin H 


z bin L 


z bin M z bin H 
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>5 


>5 >5 


>5 
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3 1 


3 
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Cluster model A y (10 -5 ) 

# (J a from fiducial, pessimistic # a a from fiducial optimistic 



z bin L z bin M z bin H z bin L z bin M z bin H 

if max = if 200 (Model 7) 5.66 2.9 4.4 1.9 3.2 4.7 2.2 

if max = l.SXioo (Model 8) 3.84 5 >6 3.4 5.6 >6 4 

Waccr = 0.28 (Model 2) 5.41 1.9 3 1.3 2.1 3.1 1.5 

oj accr = 3.5 (Model 3) 4.33 2.7 4 1.8 2.9 4.3 2 

K ph = 5 x 10 33 (Model 4) 0.70 >6 >6 >6 >6 >6 >6 

Sl = s 2 = 1.5 (Model 6) 9.39 >6 >6 >6 >6 >6 >6 



TABLE 9 

S ample models for the central Compton parameter that fall within 2-6 a from the fiducial in amplitude, described by 

IREID fc SPERGElj 120051) . THEIR BEST-FIT AMPLITUDES Ay (WHERE J/0 = Ay{M/M ) c ' M v ), and the number of a BETWEEN their 
AMPLITUDES AND THE FIDUCIAL AMPLITUDE IN THE REDSHIFT BINS AND IN BOTH THE OPTIMISTIC AND PESSIMISTIC CASES. 



Cen, R. 1997, ApJ, 485, 39 

Clowe, D., Gonzalez, A., & Markevitch, M. 2004, Astrophys. J., 
604, 596 

Clowe, D., et al. 2005, astro-ph/0511746 

da Silva, A. C, Kay, S. ±., middle, A. R,., & Thomas, P. A. 2004, 

Mon. Not. Roy. Astron. Soc., 348, 1401 
Dodelson, S. 2004, Phys. Rev., D70. 023008 
Dos Santos, S., & Dore, O. 2001, astro-ph/0106456 
Eisenstein, D. J., k. Hu, W. 1997, ftstronnvs. .1 511. 5 
Francis, M. R., Bean, R., & Kosowsky, A. 2005, astro-ph/0511161 
Hoekstra, H. 2001, A&A, 370, 743 

Hockstra, H. 2003, Mon. Not. Roy. Astron. Soc, 339, 1155 

Holder, C, Haiman, Z., & Mohr, J. 2001, Astrophys. J., 560, Llll 

Kaiser, N., & Squires, G. 1993, Astrophys. J., 404, 441 

Lapi, A., Cavaliere, A., & Menci, N. 2005, Astrophys. J., 619, 60 

Majumdar, S., & Mohr, J. J. 2004, ApJ, 613, 41 

Marian, L., & Bernstein, G. 2005, in preparation 

McCarthy, I. G., Babul, A., Holder, G. P., & Balogh, M. L. 2003, 

Astrophys. J., 591, 515 
McCarthy, I. G., Babul, A., Holder, G. P., & Balogh, M. L. 2003, 

ApJ, 591, 515 

Metzler, C. A., White, M. J., Norman, M., & Loken, C. 1999, 

Astrophys. J., 520, L9 
Motl, P. M., Hallman, E. J., Burns, J. O., & Norman, M. L. 2005, 

Astrophys. J., 623, L63 



Nagai, D. 2005, astro-ph/0512208 

Navarro, J. F., frenK, C. a., & white, S. D. M. 1997, Astrophys. 
J. 490 493 

Oh, S. P.' & Benson, A. J. 2003, Mon. Not. Roy. Astron. Soc, 342, 
664 

Ostriker, J. P., Bode, P., & Babul, A. 2005, Astrophys. J., 634, 964 
Reid, B. A., & Spergel, D. N. 2005, in preparation 
Rovchowdhurv. S.. Ruszkowski, M., & Nath, B. B. 2005, 

astro-ph/0508120 
Sheldon, a. S., et al. 2001, Astrophys. J., 554, 881 
Sheth, R. K., & Tormen, G. 1999, Mon. Not. Roy. Astron. Soc, 

308, 119 

Spergel, D. N., Verde, L., Peiris, H. V., Komatsu, E., Nolta, M. R., 
Bennett, C. L., Halpern, M., Hinshaw, G., Jarosik, N., Kogut, 
A., Limon, M., Meyer, S. S., Page, L., Tucker, G. S., Weiland, 
J. L., Wollack, E., & Wright, E. L. 2003, APJS, 148, 175 
Sunyaev, R. A., & Zeldovich, Y. B. 1980, ARAA, 18, 537 
White, M. J., Hernquist, L., & Springel, V. 2002, Astrophys. J., 
579, 16 

Wright, C. O., & Brainerd, T. G. 2000, ApJ, 534, 34 



